Detectability of a phantom-like braneworld model 
with the integrated Sachs- Wolfe effect 
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We study a braneworld model in which a phantom-like behaviour occurs with only cold dark 
matter and a cosmological constant, due to a large distance modification of gravity. With the 
addition of curvature, the geometrical tests are not strict enough to rule out models in which 
gravity is modified significantly on large scales. We show that this degeneracy in the parameter 
space is broken by the structure formation tests, such as the integrated Sachs- Wolfe effect, which 
can probe general relativity on large scales. 
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I. INTRODUCTION 

The late time acceleration of the Universe is one of 
the biggest problems in cosmology. In the framework 
of conventional general relativity, the expansion of the 
Universe at late times is dominated by a dark energy with 
negative pressure and equation of state w = pj p < —1/3. 
Several current observations suggest w < — 1, which is 
often called phantom dark energy, although the fiducial 
LCDM model with w = — 1 is still preferred if we combine 
all the data sets [HQ. From a theoretical point of view, it 
is extremely difficult to realise dark energy models with 
w < — 1: the easiest way to obtain such a model is to 
consider a ghost scalar field with the wrong sign for the 
kinetic term, although this leads to the instability of the 
vacuum There are a few successful models that lead 
to w < — 1 without having theoretical pathologies [J, Q ; 
among them, we focus on a braneworld model proposed 
by Sahni and Shtanov [6| and further developed by Lue 
and Starkman Q- 

This model is based on the Dvali-Gabadadze-Porrati 
(DGP) model of the 5D braneworld where we are sup- 
posed to live on a 4D brane in the 5D Minkowski space- 
time Q . The 4D gravity on the brane is recovered by the 
induced 4D Einstein-Hilbert action on the brane. In this 
model there are two branches of the solutions @ : in the 
first branch, known as self-accelerating, the late time ac- 
celeration can be realised without introducing any dark 
energy, while in the other, known as the normal branch, 
a cosmological constant is needed to explain the late time 
accelerated expansion of the Universe; nevertheless, the 
extra-dimensional effects modify gravity on large scales 
and the model deviates from the standard LCDM. In par- 
ticular, at the background level, the Universe behaves as 
if there were a phantom-like dark energy w < — 1. 

Besides the fact that this model mimics a phantom 
behaviour, it is known to be free of ghosts and thus rep- 



resents a healthy modified gravity theory. This is in con- 
trast with the self-accelerating branch of the DGP model 
(hereafter sDGP) where there exists a ghost at the lin- 
earised level (for a review see [l(|). Another advantage 
of the model is that there is a mechanism to recover gen- 
eral relativity on small scales. Thus with this model we 
can modify gravity on large scales significantly without 
spoiling the success of general relativity on the solar sys- 
tem scales, providing the basis for the test of the large 
distance modification of general relativity. 

In this paper, we study the phenomenological con- 
sequences of the normal branch DGP model (hereafter 
nDGP). We first present in Section II the geometrical 
tests on nDGP, looking for a parameter space which can 
be tested from structure formation, which is summarised 
in Section III. Then we present the ISW-galaxy correla- 
tions as a powerful tool to distinguish between LCDM 
and nDGP models in Section IV. Section V is devoted to 
the conclusion. 



II. GEOMETRICAL TESTS 

The cosmic expansion of the nDGP model depends on 
the usual 4D FRW metric plus the gravitational effect 
of the 5D bulk on the brane. The cosmic acceleration 
is then introduced by the brane tension, which works 
as a cosmological constant on the brane. The gravity 
at large scales is modified by the 5D gravity effects on 
the brane, which are parameterised by a transition scale 
from 4D gravity to 5D gravity. The crossover distance r c 
is defined as the ratio between 4D and 5D Planck mass 
scales 
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where M4 and M5 are the Planck scales in the 4D and 5D 
spacetime respectively. The late time expansion history 
is determined by two free parameters, the cosmological 
constant (or brane tension) A and the crossover distance 
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The Friedmann equation for an nDGP model with cur- 
vature K = — OfciJQ is given by 
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and the dimensionless expansion history E(a) is defined 
by 
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where fl\ = A/3Hq and f2 rc = {^H^r 2 ) 1 , which satis- 
fies 



&2 *ti ~T~ i2 / 



(4) 



The free parameter r c can range in theory from to the 
infinity; however, it has been shown that the deviations 
from general relativity on solar system scales are also 
controlled by r c , and the current constraints require that 
r c > Hq 1 . We can see that if r c approaches the infinity, 
then Eq. @ converges to GR, while if r c approaches Hq , 
then the 5D gravitational effect on the expansion history 
becomes maximal. 

The modification of gravity at late time screens the 
cosmological constant and makes the effective equation 
of state less than —1. We define the effective energy 
density of dark energy p e g as [H| 
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It is clearly seen that the 5D effects make the effective 
dark energy density p e g smaller. From the continuity 
equation of p e g 



p c s + 3H(1 + w e g)p e g = , 
we can derive w e s as 



(6) 
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At the current time, the effective equation of state be- 
comes 



w ef f(a = !) = -!- 



(i-n m )(n m + n A + i) : 



(8) 



where we neglected the curvature for simplicity. Provided 
that Q m < 1, we have the phantom behaviour w c g < — 1. 
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FIG. 1: Geometrical test on the nDGP by using SN+CMB+// 
observations. There is correlation observed in the projected £7^ 
and f^A plane after marginalisation of all other cosmological pa- 
rameters. 



We revisit the geometrical test on the nDGP [111 03 • 
The geometrical test on the nDGP with a flat curva- 
ture prior is not in favour of the cases for the signifi- 
cant screening effect, which rules out observable modi- 
fied gravity effects in the nDGP. However we find that 
measurable screening effects are allowed with the inclu- 
sion of curvature. We exploit the leverage arm in the 
geometrical tests at both ends of low and high redshifts. 
At low redshifts, we use the Gold SN data set [l3(. At 
high redshifts, we fix the distance to the last scattering 
surface at z\ ss = 1088^2 by fitting the harmonic space 
scale of the acoustic peak l A = i02t_^\ and matter den 



+0.0072 
-0.0095 



[HI ]. In addition to that, we 



sity fl m h 2 = 0.1268 
constrain the expansion constant Hp with the Hubble 
constant measurement, Hq = 72^g [l5j . 

With a fixed CMB prior of Q m h 2 , best fit values for w 
and Ho are correlated with each other. The theoretical 
models predicting w < — 1 have a smaller best fit value 
for Hq compared with LCDM (w = —1). Since the mea- 
sured comoving distance to zi ss is consistent with a best 
fit value for Ho in flat LCDM, the comoving distance to 
z\ ss in phantom-like braneworld models becomes longer 
than the measured distance. This worse fit for the large 
distance measured by CMB in the models with w < — 1 
can be cured by introducing a positive curvature which 
makes the distance shorter without significantly affecting 
the fit for the shorter distance measured by SNe. Con- 
sequently, a larger f^A, which realises larger screening 
effects and w < —1, is allowed with a positive curvature 
(Ofc < 0) as is shown in Fig. [TJ Hence if the curvature 
is added, there appears a degeneracy in the geometrical 
tests and the models with large modified gravity effects 
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are allowed. This degeneracy can be broken by the struc- 
ture formation test. 



III. STRUCTURE FORMATION TESTS 

There are three regimes of gravity in the nDGP model 
on different scales. On super-horizon scales, gravity is 
significantly influenced by 5D effects. In this regime, we 
cannot ignore the time evolution of metric perturbations 
and the dynamical solutions should be obtained by solv- 
ing the 5D equations of motion. The dynamical solutions 
have been obtained in the following two methods in the 
literature: a first derivation is obtained by the scaling 
ansatz in the sDGP [l6| and in the nDGP [Tjj, and the 
other is found from the full 5D numerical simulations 
[li ]. It has been shown that both approaches give iden- 
tical results, and the solutions for the perturbations are 
shown to be insensitive to the initial conditions for the 
5D metric perturbations. 

On sub-horizon scales, we can ignore the time depen- 
dence of the metric perturbations and the quasi-static 
approximations can be used [l^, l2Cj . Even on scales 
smaller than r c , gravity is not described by general rela- 
tivity due to an extra scalar degree of freedom introduced 
by the modification of gravity. In this regime, gravity can 
be described by a Brans-Dicke theory and the growth of 
structure becomes scale independent. 

We use the Newtonian gauge 



ds 2 



(1 + 2*) dt 2 + a(t) 2 (l + 2$) Sijdx^x 3 , (9) 



to describe the metric perturbations. Fig. [2] shows the 
behaviour of metric perturbations <&_ = (<£> — v I / )/2 which 
determines the integrated Sachs- Wolfe (ISW) effect both 
for the dynamical solutions and scaling solutions, for the 
models of Table[IJ In the literature, the spatial curvature 
was not introduced in the calculations, and thus we derive 
the quasi-static solutions with curvature in Appendix A. 
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TABLE I: Details of the cosmological models used. The other 
parameters are for all models: scalar spectral index n s = 0.95, 
optical depth r = 0.10, and amplitude of the primordial scalar 
perturbations A = 2.04- 10~ 9 at a pivot scale k = 0.05 Mpc -1 . 



Finally, once the non-linearity of density perturbations 
becomes important, the theory approaches general rela- 
tivity pjj This transition to general relativity is 
crucial to satisfy the tight constraints from the solar sys- 
tem experiments [22L [23[ , and will play a crucial role for 
weak lensing measures. On the other hand, for the ISW 
effect, we can safely ignore the non-linear physics. 




FIG. 2: We plot the solutions of structure formation of three nDGP 
models in the 1 — a contour of Fig. [T] compared with the LCDM 
(dotted line). Solid curves represent the quasi-static solutions of 
nDGP models with different f^, and the dashed curve attached 
to each solid curve represents the dynamic solution of each nDGP 
model at k = 10~ 3 Mpc -1 . 



The dynamical solutions are relevant to the scales of 
the large scales CMB anisotropics. We have checked that 
the difference in the large scales CMB anisotropies from 
LCDM are small given the constraints from the geometri- 
cal tests because, due to the large cosmic variance, we do 
not expect that the CMB anisotropies on these scales can 
give strong constraints on the models. The quasi-static 
solutions are relevant to the scales of ISW-galaxy cross- 
correlations. In the next section, we will study how they 
can be used to break the degeneracy that arises from the 
geometrical tests. 



IV. ISW-GALAXY CORRELATIONS 

The gravitational potential well is shallower in the 
nDGP model than in the LCDM model due to the modifi- 
cation of gravity. This is the opposite from what happens 
in the self-accelerating models [24[ where the gravita- 
tional potential well is deeper than in LCDM. The nDGP 
model predicts an earlier variation of the gravitational 
potential than the LCDM model. By cross-correlating 
galaxies at different redshifts with the CMB, one can in 
principle trace the redshift history of the decay of the 
potential. Furthermore, the cross-correlation arises from 
the well understood quasi-static (QS) regime of nDGP 
(solid curves in Fig. 

The cross-power spectrum of the CMB and a set of 
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FIG. 3: Measurement of the cross-correlation functions between six different galaxy data sets and the CMB, reproduced from The 
curves show the theoretical predictions for the ISW-galaxy correlations at each redshift for the LCDM model (black, dashed) and the three 
nDGP models of Fig. [2] (green, solid), which describe the 1 — a region of the geometry test from Fig. [l] 



galaxies <7, is given by [HI, H(| 

c;^ i ,Jfi; mrm k ' p ^ (M) . do) 

where (fc, 0) is the power spectrum of at the 

present time and the kernel ij is 
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The window function is given by 
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where the galaxy kernel I 9i is 
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(11) 



(12) 



(13) 



Under the QS approximation, the window function be- 
comes 



W 9i {k,z) = 



2 k 2 ni{z)bi(z) $_(fc,z) 

3a^#f i + z $_(fc,o) 



(14) 



where rii(z) is the redshift distribution of the galaxies 
normalised to J dzrii = 1 and bi(z) is the galaxy bias. 

First, we investigate the current status of the observa- 
tions using the data set obtained in ||, which consists 
in the measurement of the ISW effect via the real space 
cross-correlation functions (CCF) between six different 
galaxy catalogues and the CMB. The redshift distribu- 
tions of the different catalogues are partly overlapping, 
but each data set i is characterised by a median redshift 
0.1 < z~i < 1.5 around which each particular measured 
CCF is getting the biggest contribution. Thus, this anal- 
ysis represents a first step towards an exploration of the 
redshift evolution of the potentials and, ultimately, of 
gravity itself as described by (27| . 

We reproduce in Fig. [3] the measured CCF for the six 
galaxy catalogues from [2|, in order of increasing redshift: 
2MASS (excluding the small scale contaminated data), 
the main galaxy sample from the SDSS, the SDSS Lumi- 
nous Red Galaxies, NVSS, HEAO and the SDSS quasars, 
with the relative error bars which should be remembered 
are highly correlated. Looking at the theoretical curves 
in Fig. [3l we can see that the nDGP models have a very 
different prediction from the LCDM for the CCF at high 
redshift. This is in agreement with their peculiar po- 
tential evolution: the rise in the potential at high 
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FIG. 4: The theoretical predictions for the ISW-galaxy cross power spectra at each redshift for the LCDM model (black, dashed) and the 
three nDGP models of Fig. [2] (green, solid), which describe the 1 — a region of the geometry test from Fig. [T] 



redshift produces an expected negative CCF, while the 
following steeper decay leads to a positive CCF which 
becomes eventually higher than the LCDM one. 

However, it is clear that these predictions represent a 
poor fit to the high redshift data. Remembering that 
all three nDGP models in Fig. [3] are inside the la region 
from the geometry test of Section II, we can qualitatively 
see that the ISW test will produce stricter constraints by 
noticing e.g. that the quasar CCF alone has a significance 
level of 2(7, which means that at least two of the nDGP 
models will be excluded at above this level. 

Then, we study the best possible constraints which can 
be obtained by this technique with future surveys. For 
definiteness, we assume that the galaxy sets come from 
a net galaxy distributions of 

n g (z) oc zV (2/1 ' 5)2 , (15) 

where the normalisation is given by the LSST expectation 
of 35 galaxies per arcmin 2 . For the subsets of galaxies, we 
assume that this total distribution is separated by photo- 
metric redshifts which have a Gaussian error distribution 
with rms a(z) = 0.03(1 + z). The redshift distributions 
are then given by [2jj 

n i( z ) — -^-n a (z) erfc ( ^-J- | — crfc ( — — ] , 

where erfc is the complementary error function and Ai is 
determined by the normalisation constraint. 

We show in Fig. 0] the predicted cross power spectra 
obtained using this redshift tomography for the models 
of Table HI The theoretical possibility to distinguish be- 



tween them is given by the signal to noise ratio 

=Vfsky(2J+l) [C ' g]2 . , (16) 

where Cj T is the temperature power spectrum. This is 
summarised in Table HT1 
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2.6 
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0.96 


0.22 


1.5 


2.9 



TABLE II: Theoretical signal to noise ratio for the models of 
Table U with f sky = 1. 



Although the geometrical test is not able to easily 
break the degeneracy between curvature and the screen- 
ing effect, the alternative consequence for the structure 
formation by the screening effect is measurable from the 
ISW-galaxy cross-correlations. The screening of the cos- 
mological constant in nDGP2 and nDGP3 becomes effec- 
tive before the decay of the growth factor which occurs 
when the matter component becomes subdominant. This 
early screening enhances the growth factor which makes 
the potential grow. This generates anti-correlations 
in the ISW-galaxy cross-correlations at high redshifts, 
which leaves observable signatures as is shown in Fig. 2] 
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From Table HT1 it is expected that this effect on the struc- 
ture formation can be observed at around 50% noise level 
for nDGP2 and 25% noise level for nDGP3. This is an il- 
lustration how we can break the degeneracy between cur- 
vature and the screening effect in the geometrical tests 
by using the structure formation tests. 

V. CONCLUSION 

In this paper, we studied the observational constraints 
on the normal branch DGP model in which a phantom- 
like behaviour occurs only with cold dark matter and a 
cosmological constant. The geometrical tests using the 
gold SN data set, CMB and the HST key project are 
not enough to rule out models in which gravity is signifi- 
cantly modified on cosmological scales. We then showed 
that the structure formation tests performed using the 
integrated Sachs- Wolfe (ISW) effect can break the de- 
generacy in the parameter space. 

The current measurements of the ISW effect obtained 
in @ are indeed as competitive as the geometrical tests. 
This is due to the fact that, in the nDGP model, the 
cross-correlation with galaxies becomes negative at high 
redshift due to the peculiar behaviour of the metric per- 
turbations caused by the modification of gravity. This 
demonstrates that the structure formation tests are very 
promising tools to distinguish between general relativity 
and modified gravity models. We also showed that it is 
possible to track the evolution of the potentials by cross- 
correlating the ISW with galaxies at each redshift in fu- 
ture observations. It is very likely that in the future the 
ISW effect will provide one of the strongest constraints 
on the model. We will present the full likelihood analysis 
using the latest data sets in a forthcoming paper. 
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APPENDIX A: QUASI-STATIC SOLUTIONS 
WITH CURVATURE 



In the Gaussian normal coordinates, the 5D metric is 
given by [13] 

ds 2 = dy 2 - n(y,t) 2 dt 2 + a(y,t) 2 S lj dx l dx :i , (Al) 
where 



a{y,t) = a(t) 
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The extrinsic curvature of the brane is determined by 
the first derivative of the metric with respect to y at the 
brane (y = 0): 
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Defining the comoving density perturbations 
pA = 8p — 'SHaSq, 
the Poisson equation is obtained as 
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The traceless part of the space-space component of the 
effective Einstein equations gives 
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r c [{a 1 /a) + (n'/n)] 
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(A8) 



The Weyl density perturbations SpE, 5qE and Site should 
be determined by the constraint equations 
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The constraint equations are not closed and we need 
additional information by solving the 5D equation of mo- 
tion. In the quasi-static limit, we can impose the condi- 
tion on SpE and Stye from the bulk equation as [20] | 



5pE = 2k 2 SnE- 
Then the constraint equations give 

-l + {a'/a)r c + (n'/n)r c 



SpE 



-3 + 4(a'/a)r c + 2{n'/n)r c 



(All) 



pA, (A12) 



and SqE = 0. The Poisson equation and the traceless 
part of Einstein equations give 
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where 



m = i - 



which can be written as 



/3(t) = 1+ 2H 2 r c H 2 - 
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